Concentration of solutions for a 
singularly perturbed mixed problem 
in non smooth domains 



Serena Dipierro 



SISSA, Sector of Mathematical Analysis 
Via Bonomea 265, 34136 Trieste, Italy 
E-mail address: dipierro@sissa.it 



Abstract 



We consider a singularly perturbed problem with mixed Dirichlet and Neumann boundary 
conditions in a bounded domain C R" whose boundary has an (n — 2)-dimensional sin- 
gularity. Assuming 1 < p < ^J^f, we prove that, under suitable geometric conditions on the 
boundary of the domain, there exist solutions which approach the intersection of the Neumann 
and the Dirichlet parts as the singular perturbation parameter tends to zero. 
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1 Introduction 

In this paper we study the following singular perturbation problem with mixed Dirichlet and 
Neumann boundary conditions in a bounded domain 51 C M" whose boundary dfl is non smooth: 



Here p G {^1, ^Jzfj is subcritical, v denotes the outer unit normal at 90 and e > is a small 

parameter. Moreover d^^, do^ are two subsets of the boundary of such that the union of 
their closures coincides with the whole 9f2, and their intersection is an (n — 2)-dimensional smooth 
singularity. 

Problem ([T]) or some of its variants arise in several physical or biological models. Consider for 
instance the study of the population dynamics: suppose that a species lives in a bounded region 
whose boundary has two parts, On^, which is an obstacle that blocks the pass across, and Od^, 
which is a killing zone for the population. Moreover ([T]) is a model of the heat conduction for 
small conductivity, when there is a nonlinear source in the interior of the domain, with combined 
isothermal and isolated regions at the boundary. 

Concerning reaction- diffusion systems, this phenomenon is related to the so-called Turing's 
instability. More precisely, for single equation with Neumann boundary conditions it is known 




e^Au + u = uP in Q, 
- = on On^, u = on Od^, 
% in VI. 



(1) 




that scalar reaction-diffusion equations in a convex domain admit only constant stable steady state 
solutions; see |CHj . [Mat| . On the other hand, as noticed in (Tuj, reaction-diffusion systems with 
different diffusivities might generate non-homogeneous stable steady states. A well-known example 
is the Gierer-Meinhardt system, introduced in |GM| to describe some biological experiment. We 
refer to [Nlj, |NTYj for more details. 

Another motivation comes from the Nonlinear Schrodinger Equation 

where h is the Planck constant, V is the potential, and 7 and m are positive constants. In fact, 
if we analyze standing waves and consider the semiclassical limit fi. — > 0, we obtain a singularly 
perturbed equation; see for example |ABC| . |AMj . |AMSj . [FWj . and references therein. 

Let us now describe some results which concern singularly perturbed problems with Neumann 
or Dirichlet boundary conditions, and specifically 

— e^Au + u — u'P in fi, 

= on an, (2) 



u > in n 



and 




in ri, 

on dfl, (3) 
in n. 

The study of the concentration phenomena at points for smooth domains is very rich and has 
been intensively developed in recent years. The search for such condensing solutions is essentially 
carried out by two methods. The first approach is variational and uses tools of the critical point 
theory or topological methods. A second way is to reduce the problem to a finite-dimensional one 
by means of Lyapunov-Schmidt reduction. 

The typical concentration behavior of solution Uq,^ to ([T]) is via a scaling of the variables in 
the form 

'x - 



UQ^,{x)^Ui^^^j, (4) 

where Q is some point of H, and C/ is a solution of the problem 

-AU + U = UP in M" (or in M'| = {{xi, . . . , x„) e M" : a;„ > 0}), (5) 

the domain depending on whether Q lies in the interior of fl or at the boundary. When p < 
(and indeed only if this inequality is satisfied), problem ([5| admits positive radial solutions which 
decay to zero at infinity; see |BLj . [St| . Solutions of ([I]) with this profile are called spike-layers, 
since they are highly concentrated near some point of il. 

Consider first the problem with pure Neumann boundary conditions. Solutions of ^ with 
a concentration at one or more points of the boundary dCl as e — )■ are called boundary- spike 
layers. They are peaked near critical points of the mean curvature. In particular, it was shown 
in [NT1| . |NT2| that mountain-pass solutions of ^ concentrate at 9ri near global maxima of the 
mean curvature. One can see this fact considering the variational structure of the problem. In 
fact, solutions of ([2| can be found as critical points of the following Euler-Lagrange functional 



{u) = l f {e^\Vu\^ + u^) dx ^— f \u\P+^dx, ueH^in) 
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Plugging into I^^n a function of the form Q with Q E dil one sees that 

I.,N {Uq,,) = Coe" - Cie"+'H (Q) + o (e"+i) , (6) 

where Cq, Ci are positive constants depending only on n and p, and H is the mean curvature; see 
for instance |AMj . Lemma 9.7. To obtain this expansion one can use the radial symmetry of U and 
parametrize 9f2 as a normal graph near Q. From the above formula one can see that the bigger 
is the mean curvature the lower is the energy of this function: roughly speaking, boundary spike 
layers would tend to move along the gradient of H in order to minimize their energy. Moreover one 
can say that the energy of spike-layers is of order e", which is proportional to the volume of their 
support, heuristically identified with a ball of radius e centered at the peak. There is an extensive 
literature regarding the search of more general solutions of ^ concentrating at critical points of 
H; see [DFW] . [G?] . pFW] . [G^, [Li], pJTj . [NPT] . |We) . 

Consider now the problem with pure Dirichlet boundary conditions. In this case spike layers 
with minimal energy concentrate at the interior of the domain, at points which maximize the 
distance from the boundary; see |LNj . |NW| . The intuitive reason for this is that, if Q is in the 
interior of fl and if we want to adapt a function like Q to the Dirichlet conditions, the adjustment 
needs an energy which increases as Q becomes closer and closer to dil. Following the above 
heuristic argument, we could say that spike layers are repelled from the regions where Dirichlet 
conditions are imposed. 

Concerning mixed problem ([I]), in two recent papers [GMMPlj . |GMMP2j it was proved that, 
under suitable geometric conditions on the boundary of a smooth domain, there exist solutions 
which approach the intersection of the Neumann and the Dirichlet parts as the singular perturba- 
tion parameter tends to zero. In fact, denoting by u^^q an approximate solution peaked at Q and 
by de the distance of Q from the interface between and do^, then its energy turns out to be 
the following 

/, (wQ,e) = Coe" - C^e^+^H (Q) + e"e"2^(i+°(i» + o {e"+^) , (7) 

where 1^ is the functional associated to the mixed problem. Note that the first two terms in ([t]) are 
as in the expansion ([6]), while the third one represents a sort of potential energy which decreases 
with the distance of Q from the interface, consistently with the repulsive effect which was described 
before for 

In almost all the papers mentioned above the case of f2 smooth was considered. Concerning 
instead the case of fl non smooth, in [Dij the author studied the concentration of solutions of 
the Neumann problem ^ at suitable points of the boundary of a non-smooth domain. Assuming 
for simplicity that 51 C K'^ is a piecewise smooth bounded domain whose boundary d^l has a 
finite number of smooth edges, one can fix an edge F on the boundary and consider the function 
a : F — )■ K which associates to every Q GT the opening angle at Q, a (Q). Then it was proved that 
this function plays a similar role as the mean curvature H for a smooth domain. In fact, plugging 
into J^^iv a function of the form (jll) with Q £ T, one obtains the analogous expression to (j6| for 
this kind of domains, with Coa{Q) instead of Co. Again, one can give an heuristic explanation 
considering the fact that in this case one has to intersect the ball of radius e, which is identified 
with the support of the solution, with the domain, obtaining the dependence on the angle a {Q). 

We are interested here in finding boundary spike-layers for the mixed problem ([I]). We call F 
the intersection of the closures of djqVL and d^fl, and suppose that it is an (n — 2)-dimensional 
smooth singularity. Moreover we denote by H the mean curvature of dfl restricted to the closure 
of On^, that is H : d^i^ M. 

The main result of this paper is the following: 
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Theorem 1.1. Let fl C M", n > 2, be a bounded domain whose boundary has an {n — 2)- 

dimensional smooth singularity, and 1 < p < ^j^l (1 < p < +oo if n = 2). Suppose that diy^l,, 

are disjoint open sets of dfl such that the union of the closures is the whole boundary of 
and such that their intersection T is the singularity. Suppose Q d T is such that a (Q) ^ and 
7J|r is critical and non degenerate at Q, and that V H (Q) ^ points toward Od^- Then for e > 
sufficiently small problem ([ij admits a solution concentrating at Q. 

Remark 1.2. (a) The non degeneracy condition in Theorem \l.l\ can be replaced by the condition 
that Q is a strict local maximum or minimum of H\r, or by the fact that there exists an open 
set VofT containing Q such that H (Q) > snpgy H or H (Q) < infgy H. 

(b) With more precision, as e 0, the above solution possesses a unique global maximum point 
Qe & Qn^, and dist {Q^,r) is of order e log ^. 

The general strategy for proving Theorem |1.1| rehes on a finite-dimensional reduction; see for 
example the book [AMj . One finds first a manifold Z of approximate solutions to the given problem, 
which in our case are of the form Q , and solve the equation up to a vector parallel to the tangent 
plane of this manifold. To do thi s on e can use the spectral properties of the linearization of ([s]), see 
Lemma 



4.3 Then, see Theorem ] 2. 2 [ one generates a new manifold Z close to Z which represents 
a natural constraint for the Euler functional of Q, which is 

H - ^ / (e'lVup + u') dx - f \u\P+'dx, ueHh m , 

where Hj-, (51) is the space of functions (51) which have zero trace on Sofl. By natural constraint 
we mean a set for which constrained critical points of are true critical points. 

Now, we want to have a good control of the functional \^. Improving the accuracy of the 
functions in the original manifold Z, we make Z closer to Z \ in this way th e main term in the 



constrained functional will be given by 1^ |z, see Propositions 3.12 3.14 3.15 To find sufficiently 
good approximate solutions we start with those constructed in literature for the Neumann problem 
([2| (see Subsection 2.2) which reveal the role of the mean curvature. The problem is that these 



functions are non zero on Sjjfl, and even if one use cut-off functions to annihilate them the 
corresponding error turns out to be too large. Following the line of |GMMPl] and |NW| . we will 
use the projection operator in _ff^ (H), which associates to every function in this space its closest 
element in H]-, (fl). To study the asymptotic behavior of this projection we will use the limit 
behavior of the solution L/ to ([5|: 

lim e^'r^U (r) = c„,p, (8) 

r— >-+oo 

where r = and c„^p is a positive constant depending only on the dimension n and p, together 
with 

U' ir ) U" (r) 

lim 7^ = - lim ^i-i^^-l, (9) 

r->+oo U [r) r^+oo (J (r) 

as it was done in some previous works, see for instance [LN| and |We2] . Moreover, we will work 



at a scale d ~ e|loge|, which is the order of the distance of the peak from F, see Remark 1.2 
(b). At this scale both and du^ look fiat; so we can identify them with the hypersurfaces 
of equations a;„ = and Xi tana + Xn = 0, and their intersection with the set {xi = x„ = 0}. 
Note that a = a (Q) is the angle between xi and Xn at a fixed point Q G T. Then we can replace 
fl with a suitable domain E^i, which in particular for < a < tt is even with respect to the 



coordinate Xn, see the beginning of Subsections 3.1 and 3.2 Now, studying the projections in this 



domain, we will find functions which have zero x„-derivative on {xn — 0}\9S£), which mimics the 
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Neumann boundary condition on On^- After analyzing carefully the projection in Subsections 3.1 
3.2 we will be able to define a family of suitable approximate solutions to ([T]) which have sufficient 
accuracy for our analysis, estimated in Propositions |3. 12} |3.14[ [3.15[ 

We can finally apply the above mentioned perturbation method to reduce the problem to a 
finite-dimensional one, and study the functional constrained on Z. We obtain an expansion of the 
energy of the approximate solutions, which turns out to be 



/, (ue.g) = Coe" - (7ie"+iff (Q) + e^e'^- 
in the case < a < ^, and 

/, K,q) = Coe" - Cie^+'H (Q) + e"e-2^(i+°(i)) + o {e^+') , 

in the case f < a < 27r. As for ([?]), we have that the first two terms come from the Neumann 
condition, while the others are related to the repulsive effect due to the Dirichlet condition. Let us 
notice that, in the first case, in the terms related to the Dirichlet condition appears the opening 
angle a, whereas in the second case it does not; this phenomenon comes from the fact that the 
distance of the point Q from the Dirichlet part d]jil depends on a only if < a < |. 

Concerning the regularity of the solution, following the ideas in |Grij , it is possible to say that 
it is influenced by the presence of the angle. In fact, the solution is at least in the interior of 
the domain, far from the angle; whereas, near the angle, one can split the solution into a regular 
part and a singular one, whose regularity depends on the value of a. For more details about the 
regularity of solutions in non-smooth domains we refer the reader to the book [Gri . 

The fact that the solution u is in the interior of the domain allows to say also that it is 
strictly positive, by using the strong Maximum Principle. In fact, we have that u > in the 
domain. Moreover, if there exists a point xq in the interior of the domain such that u (xq) — 0, we 
can consider a ball centered at Xq of small radius suet that it is contained in the domain; since in 
the ball u is we can conclude that u cannot be zero in xq. 

The plan of the paper is the following. In Section [2] we collect some preliminary material: 
we recall the abstract variational perturbative scheme and some known results concerning the 
Neumann problem ([2]). In Section [s] we construct a model domain to deal with the interface, 
analyze the asymptotics of projections in and then construct approximate solution to ([l]). 
Finally in Section |4] we expand the functional on the natural constraint, prove the existence of 
critical points and deduce Theorem 



Notation 

Generic fixed constant will be denoted by C, and will be allowed to vary within a single line or 
formula. The symbol o (t) will denote quantities for which tends to zero as the argument t 
goes to zero or to infinity. We will often use the notation d(l-|-o(l)), where o(l) stands for a 
quantity which tends to zero as d — > +oo. 



2 Preliminaries 

We want to find solutions to ([ij with a specific asymptotic profile, so it is convenient to make the 
change of variables x i— )■ ex, and study ([I]) in the dilated domain 

:= -n. 

e 
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Then the problem becomes 

—Am + u = in Hj, 

f^==0 ondN^,, u = ondo^e, (10) 

u > in rjj, 

where On^c and dn^t stand for the dilations of and respectively. Moreover we denote 
by Fe the intersection of the closures of On^e and dn^e- 



Solutions of ( 10 1 can be found as critical points of the Euler-Lagrange functional 



= {\Vu\^+u^)dx —[ \u\P+'^dx, ueHlifl^). 

Here Hj^ (il^) denotes the space of functions in (fl^) with zero trace on Od^c- 

In the next subsection we introduce the abstract perturbation method which takes advantage 
of the variational structure of the problem, and allows us to reduce it to a finite dimensional one. 
We refer the reader mainly to |AM| . |Ma| and the bibliography therein for the abstract method. 
In our case we will use some small modifications of the arguments in the latter references which 
can be found in Subsection 2.1 of IGMMPlj . 

2.1 Perturbation in critical point theory 

In this subsection we recall some results about the existence of critical points for a class of func- 
tionals which are perturbative in nature. Given an Hilbert space H, which might depend on the 
perturbation parameter e, we consider manifolds embedded smoothly in H, for which 

i) there exists a smooth finite-dimensional manifold C H and C, r > such that for any 
z e Zg, the set Z^^ n Br (z) can be parametrized by a map on Bf'' whose norm is bounded 
by C. 

Moreover we are interested in functionals : _ff — > M of class C^''*' which satisfy the following 
properties: 

ii) there exists a continuous function / : (0, eg) ^ M with limc_j.o / (e) = such that ||/g(z)|| < 
/ (e) for every z G Z^; moreover [q\\\ < f (e) ||q|| for every z € and every q e T^Z^; 

iii) there exist C, 7 G (0, 1], ro > such that < C in the subset {u : dist (u, Z^) < rg}; 

iv) letting : H ^ (T^Z^)^, for every z € Z^, be the projection onto the orthogonal complement 
of TzZ^, there exists C > 0, independent of z and e, such that PzI'J{z), restricted to (TzZ^)'^, 
is invertible from (r^Z^)^ into itself, and the inverse operator satisfies {PzI'J {z))~^ < C. 

We set W — {TzZf)^ , and look for critical points of in the form u = z + w with z ^ Z^ and 
It; G W . li Pz : H W is ds, in iv), the equation I'^{z + w) = is equivalent to the following 
system 

Pzl'e {z + w) = {the auxiliary equation) , 

{Id — Pz) I'^ {z + w) = {the bifurcation equation) . 

Proposition 2.1. (See Proposition 2.1 in IGMMPlj ) Let i) — iv) hold true. Then there exists 



eg > with the following property: for all \e\ < eo '^''^d for all z ^ Z^, the auxiliary equation in ( 11 ) 
has a unique solution w = We{z) G W , which is of class with respect to z £ Z^ and such that 
\\we{z)\\ < Ci/(e) as \e\ — >■ 0, uniformly with respect to z d Z^. Moreover the derivative of w with 
respect to z, w'^ satisfies the bound 11^^(^)11 < CCif {e)'^ . 
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We shall now solve the bifurcation equation in (11). In order to do this, let us define the reduced 
functional : — >■ M by setting le{z) — I^{z + u'c(z)). 

Theorem 2.2. (Sec Proposition 2.3 in IGMMP1\I ) Suppose we are in the situation of Proposition 



2.1 and let us assume that has, for \e\ sufficiently small, a stationary point z^. Then 



z^ + w{ze) is a critical point ofl^. Furthermore, there exist c,r > such that if u is a critical point 
of le with dist Z^^c) < r, where Z^^c ^ {z £ Z^ : dist [z, dZ^) > c}, then u has to he of the form 
Ze + for some z^ £ Z^. 

2.2 Approximate solutions for ([!]) with Neumann conditions 

In this subsection we introduce some convenient coordinates which stretch the boundary and we 
recall some results from [AM, and (GMMPl; concerning approximate solutions to the Neumann 
problem. 

First of all it can be shown that near a generic point Q G F the boundary of 51 can be described 
by a coordinate system y — [yi, ...,?;„) such that 

(a) coincides with {j/„ = 0}, 

(b) d£)^l coincides with {yi tana + y^ = 0}, where a = a (Q) is the opening angle of F at Q, 

(c) the corresponding metric coefficients are given by gij — Sij + o (e). 
For further details we refer the reader to [Dij . 

Remark 2.3. (i) We stress that, in the new coordinates y, the origin parametrizes the point Q, 
and those functions decaying as \y\ — > +oo will concentrate near Q. 

(ii) It is also useful to understand how the metric coefficients gtj vary with Q. Notice that 
condition (c) says that the deviation from the Kronecker symbols is of order e, and we are 
working in a domain scaled of j ; hence a variation of order 1 of Q corresponds to a variation 
of order e in the original domain. Therefore, a variation of order 1 in Q yields a difference 
of order in gij, and precisely 

with a similar estimate for the derivatives of the inverse coefficients g^^ . For more details 
see the end of Subsection 9.2 in [AM] . 

Suppose that this coordinate system y is defined in B^^ (Q), with /zq > sufficiently small. Now, 
in this set of coordinates we choose a cut-off function x^o with the following properties 

Xtio {^) = 1 in Btta. (Q) , 

Xmo i^) - in M" \ B^ iQ) , 

IVXmoI + IV^XmoI <C in B^{Q)\B^{Q), 

and we define the approximate solution He.Q as 

Ue^Q (y) ■■= Xa.0 (ey) (Uq (y) + ^wq (y)) , (12) 

where Uq (y) = U {y — Q) and wq is a suitable function obtained in Subsection 2.2 of |GMMPl| 
by a small modifications of Lemma 9.3 in |AMj . satisfying the following estimate 

\WQ {y) I + \^WQ (y) I + l^^WQ (y) I < Co (1 + I J/I^) e"!^! , (13) 
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where Co and K are constants depending on fJ, H, n and p. 

The next result collects estimates obtained following the same arguments of Lemmas 9.4, 9.7 
and 9.8 in [XM) . 

Proposition 2.4. There exist C,K > such that for e small the following estimates hold 



1^1 (y) < ^ ^ '^'''^ - 



Mo 

4£ ' 



dUg - ^ Ce-c7 for t < lyl < 

I - A,z.,g + - uIq\ (y) < I ^ ^ ^ ^. 



where 



h,N (i2,,Q) = Co - Ciei/ (eQ) + o {e') ; — 7^,^ (u^.q) = -Cie'H' (eQ) + o (e 



An immediate consequence of this proposition is that 

\\I'e iue,Q)\\ < Ce^ for all Q e On^c such that dist {Q, T,)> ^, (14) 
where C > is some fixed constant and /iq is as before. 



3 Approximate solutions to (10) 



To construct good approximate solutions to (10 1, we will start from a family of known functions 



which constitute good approximate solutions to (10 1 when we impose pure Neumann boundary 



conditions. Since we have to take into account the effect of the Dirichlet boundary conditions, we 
will modify these functions in a convenient way. Following the line of [GMMPlj and jNW| , we will 
use the projection operator onto Hj^ i^e), which associates to every element in its closest 

point in i^e)- Explicitly, this is constructed subtracting to any given u G (ft^) the solution 
to 




in rig, 

on C^o^^e 

on disi^e 



(15) 



This solution can be found variationally by looking at the following minimum problem 



inf 

v—u on do^e 



dx 



Instead of studying ( 15 ) directly, it is convenient to modify the domain in order that the region 
of the boundary near becomes flat. We fix Q e and consider the opening angle of F,. at Q, 
a = a (Q). Since the construction of this new domain is different for < a < tt and tt < a < 27r, 
we will study separately the two cases in the following two subsections. 
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3.1 Case < a < vr 

For technical reasons we construct a domain S in the following way: we consider two hypersurfaces 
defined by the equations Xi tana + a;„ = and Xi tana — a;„ = 0, which obviously intersect at 
{xi = Xn = 0}. Then we close the domain between the two hypersurfaces with xi<OifO<a<| 
and with a;i>0if^<a<7r with a smooth surface, in such a way that the scaled domain 



(16) 



defined for a large number D, contains a sufficiently large cube. In we denote by To the 
singularity, which lies on {xi = Xn = 0}. The following figure represents a section of the domain 
in the plane xi, Xn- 



X ' 

n 

1 " IN 








X 

1 



The advantage of dealing with this set is that if we solve a Dirichlet problem in with data even 
in Xni then for suitable boundary conditions the solution in the upper part S^) H {x„ > 0} will be 



qualitatively similar to that of (15) 



Our next goal is to consider the following problem 

-A(^ + (^ = in SdD, 
(p = \J {■ - dQo) on dT,dD, 

where Qq — (—1, 0, • • • ,0). By a scaling of variables, this problem is equivalent to 

-^A(/j + (p==0 in Ejj, 



(17) 



(18) 



3.1.1 Asymptotic analysis of (18 1 



First of all we need to know if (18 1 is solvable. It follows from Lemma 3.1 in |GMMPl] : in 
fact, making a modification of some arguments in |Grij . they construct barrier functions for the 
operators A and —A + 1 at all boundary points of the set S. This guarantees, via the classical 



If we consider the function 



Perron method, the existence of a solution for the problem (18). 

— ^ \og(p, then satisfies 

iA(/.-|V0|2 + l = O inE^, 
(/. = -ilog(t/(d(--go))) ond^D- 

Using the limit behavior of the function U given by (|8| , it is easy to show the following: 



(19) 
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Lemma 3.1. For any fixed constant D > we have that 
1 



^ log [U [d (• - Qo))) ^ I • -Qol uniformly on d^o (20) 



Since Lemma 3.1 states that the boundary datum is everywhere close to the function \x — Qo\ 



it is useful to consider the following auxiliary problem 

iA(/.-|V0|2 + l = O inSc, 
(/) = |a; - Qo\ on 5E_d. 



(21) 



Lemma 3.2. Let D > 1 be a fixed constant. Then, when d oo, problem (21) has a unique 
solution (f)'^ , which is everywhere positive, and which more precisely satisfies the estimates 

. < ix) <C in S^, (22) 

V tan^ a + 1 

?/ < a < ^ , and 

l<(t)'^{x)<C in Ed, (23) 
*/ f ^ ^ I"? where C depends only on D and E. 

Proof. Applying the transformation inverse to the one at the beginning of this subsection and 
using the existence of barrier functions for the operator —A + 1, as shown in GMMPlJ, Lemma 
3.1, we get existence. Uniqueness and positivity of follows from the maximum principle. 



W 



To prove the estimates (22| and (23), we can reason as in [GMMPl] . Lemma 3.4, or in 
Lemma 4.2. In the case < a < we have that 6''' (x) = , - in En is a subsolution to (21 ) 



since dist {Qo, d'So) — ytan"a+i ' 'whereas, in the case § < a < tt, we have that dist {Qo, dY.D) — 1 
and then the subsolution is given by 01 (x) = 1. Moreover, in both the two cases, the function 
4'+ {x) = C + Xi is a supersolution for C sufficiently large. Then our claim follows. □ 

We next show some pointwise bounds on (j>'^, which in particular imply a control on the gradient 
within some region in the boundary of E^i. We obtain gradient bounds only near smooth parts of 
the boundary, away from the singularity Td. 



Lemma 3.3. Let D > 1 be as in Lemma \3.2\ Then, there exists a constant C > such that for 
any a > sufficiently small there exist S > and d^ > so large that 

\(l)'^{x) - (j)'^ {z^) I < C|a; - z^\, G dY^o, dist {z^, DTd) > a, \x - z^\ <d,d> d„. 

In the above formula denotes the point in dY,D closest to x. 

Proof. Let us first consider the case < a < ^. Let us fix cr > small and consider, for every 
Q < 5 < 5 = cr tan a, the points x G E^i of the form z + Su {z), where z € D'Ed and v {z) is the inner 
unit normal at z. Note that there is no problem in the representation of x if dist {z , DT o) > a; 
whereas if dist{z,DT d) < cr, we follow the inner normal direction given by v {z) and stop at 
x„ = if we reach this hyperplane at a distance from the boundary smaller than 6. Let us call 
this set of points x G 'Sd at distance 6 from the boundary. Note that the A^'s are all disjoint as S 
varies in [O, d\ . Now in A^ we can define the functions 

01 {x)^ \zi {x)~Qo\+MSi {x), 

02 {x) = \z2 {x)-Qo\+MS2 {x). 
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where Zi (x) , Z2 (x) are the points in 9S u closest to x with the n-th coordinate respectively positive 
and negative; Si (x), S2 {x) give the distance of x from zi (x), zi (x). If we set 

0+ (x) = min {01 (x) , 02 (x)} , 

we choose the constant M so large that 0"^ (x) > (j)'^ (x) when x G {z + 5i' (z): z£ dT.D}. The 
existence of such constant M is guaranteed by Lemma [3^ 

Next we consider a smooth function p S (K"), such that suppp C -Bi (0), and p (x) dx = 
1. Moreover we define the function 



X{x) 



=6^ {x) + 2S (x) 



for x e {z + Siy {z) : z £ dJ^oJ £ [O, S] } 



Then we construct a moUifiers 



(24) 



in such a way that the support of each px{x) depends on the point x, and, in particular, it shrinks 
to a point when we are close to the boundary. 

Finally we regularize (j)'^ using the convolution with the mollifiers defined in ( 24 1 . Then we 
obtain the following smooth function 

<l>+ (x) = * Px{-}) i^) = / ^iix-y) Px(x) (y) dy. 
It is easy to see that, for i = 1, . 



,n - 1, 



dxi 



< 0(1) 



1 



zM. 



(25) 
(26) 



Moreover, using ( 25 1 and ( 26 ) , we have that 



dxi 



dxi 



^{x-y) px(^) (y) dy+ f 0^ (x - y) ^^^^ {y) ^ (x) dy 



dx,. 



< o{l)-^M+ I J'i(x-y)^{y)^{x)dy. 



dp\(x) 



9A 

dxi 



(27) 



Now we need an estimate for the last term in (27), let us call it A. If we add and subtract 0'J_ (x) 
in the integral, we obtain 



A 



/ i'i{x-y)-4>'i{x) 



dxi 



(x) 



01 (x - y) - 0^ (x) 



dp\{x) 
dX 



(y) dy, 



in the last step we have used the fact that 0^ (x) and ^ (x) do not depend on y, and the fact 

that Jj,„^^Mr^ (y) dy = ^ px(x) (y) dy ^ 0, since px[x) (y) dy = 1, for every A > 0. Now, 
from ( |24[ ), a simple computation yields 

dp\{x) 



dX 



{y) = -nX-^^-'{x)p 



A(x) 



- A" 



(x) yVp 



A(x) 
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Then, using the fact that ^ (x) ~ —Cxi, for some positive constant C, and making the change 
of variable y = X (x) z, we have 



A = CX-^ {x) Xi / te {x-X (x) z) - 4>'l (x)] ■ [p (z) + zS/p (z)] dz. 



Since is a Lipschitz function, from (28 1 we get that 



(28) 



A<Cxi \z\-[p{z) + z\'p{z)]dz, 

and then A < o{l). It follows that, for M sufficiently large, the norm of V0^[ can be arbi- 
trarily big on its domain. By (22), if M is large then (j)'^ is everywhere bigger than on 
I]_D n {dist (•, d'Eo) = S}, so (j)'^ is a supcrsolution of (21) in S/j n {dist (•, 9Sd) < S]. 

On the other hand, we claim that the function (j>zjx) — \x — Qq\ is a subsolution of (21) in 
n {dist (-jSEu) < In fact, if we consider the set Ed \ Bg^^.^ (Qo), where 6 (d) is a small 
positive number depending on d, we can see by easy computation that here 0l satisfies 



1 = 



n-1 
d\x ~ Qo\ ' 



< 



Hence we 



Moreover, since (j>'^ is positive, we can choose S (d) sufficiently small so that 
obtain that cjy'L < (j)'^ in the closure of Ed n {dist (•, 9Ed) < 

Finally, the conclusion follows from the fact that 0l and (p'l_ coincide on the set 

{x edT,D ■■ dist (x, DTd) > cr} 

and that we have uniform bounds on the gradient here, independently on d. 

In the case ^ < a < tt, we can repeat essentially the same construction of the proof of Lemma 
3.5 in [GMMPlj and obtain the same conclusion. □ 



Using the same arguments as in Lemma 3.6 in jGMMPlj we are able to extend the gradient 
estimate which follows from the previous lemma to a subset of the interior of the domain. 



Lemma 3.4. Let D > 1 be as in Lemma \3.2\ Then, there exists a constant C > such that for 
any a > sufficiently small there exists d„ > Q so large that 



(x) I < C in {x e Sr, : dist (x, DT £)) > cr} , d > d^ 



(29) 



The next proposition is about the asymptotic behavior of the solutions of (21 ). 
Lemma 3.5. Let (p'^ be the solution of (21], then we have that 
cj>^ {x)^^{x) 



inf (Ix — z| + |z — Qol) J as d — oo. 



(30) 



uniformly on the compact sets of Yi^. 



Proof. We will show ( 30 ) in two steps: 



1) we prove that the function on the right-hand side of (30) is the supremum of all the elements 
of 

F^{ve W^'°° (Sd) ■■v{x)< \x-Qq\ on dY.uA^v\ < 1 a.c. in Sd}; 
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2) we prove that for any sequence — )■ oo, there is a subsequence dk, — >■ oo such that — >■ ^ 
uniformly on the compact sets of E^i as dk, oo. Then it follows that ^ uniformly on 

the compact sets of S/j as d ^ oo. 

We first prove 1). To begin we show that (f) G F. If xi,X2 & and Z2 & dHo realizes the infimum 
for X2, we have 

\4>{xi) -4>{X2) I < ||a;i - 02! + 1-22 - Qol - \X2 - Z2\ - \Z2 - Qoll < \XI - X2\. 

Then, taking xi,X2 close, we get </> G W^'°^ (^d) and |V0| < 1 a. c. in Sjj. Moreover, it is easy 
to see that (j) {x) = \x — Qq\ ii x €: dT,D- We next show that is the maximum element of F. We 
construct a 5 neighborhood Sfj of Sjj in this way: consider Qo = (—1,0, ••• ,0) and, for every 
z G the line from Qq to 2. If 5 > is small enough, each point x in E'^ \ E/3 is uniquely 

determined by the equation x = z -\-~Er{z), where 2; G SE^ is the intersection point of the line 
from Qo to x with dY,D, r (z) is the unit outer vector on the line, and < J < ^^^^^^^ ; here 9 (z) is 
the angle between r (z) and the unit outer normal at z, v {z), in the plane generated by r (z) and 
V [z). Note that for the point on the boundary z G {2:1 = 2;„ = 0} we can consider v (z) just taking 
the normal to the hypersurface defined by the equation x\ tan a + a;„ = or to the one defined by 
the equation xi tana — x„ = 0, and it is well defined since the angle 6 {z) is the same for those 
points. In addition, the map x — > (z, (5) is continuous in E^ \ E^j. 

Now, we can extend every w G F to a -u G W^'°^ (Ef,), taking V = V Yjz) and v {x) = v {z) for 
X G E^ \ E/j. Moreover, if we consider the function 



K{x) 



1 in Eb, 

1 + C5 inEi,\ED, 



for some large constant C > independent of 5, we get jV^I < if a. e. in E^. Now, we regularize 
V using the convolution with mollifiers, that is considering, for A > small enough, v\ := v * px, 
with px {x) = A-"-p (x/A), pGC^ (M"), suppp c Bi (0), J^^ p (x) dx = 1. Then we have 

IVi^aI < |V^| *px<K*Px<^ + CX 

on Ed and vx ^ ?; in C (E^j) as A — >■ 0. Let now x,y & Ejj and consider the function ^ {t) = 
tx + {1 — t) y, for t G [0, 1]; then we can estimate 

I^^A {x) - VX (y) I < \Vvx (*)) I • \^\dt < J\'^ + CM ' k - y\dt < (1 + CX) ■ \x - y\. 

Letting A — >■ 0, we obtain |?; (x) — v{y)\ < \x — y\. Hence v{x) < v (y) + \x — y\, and v (x) < 
\y — Qa\ + \x — y\ for all y G dT,]j. So i; < 0. 

We next prove 2). By gradient estimate and the Ascoli-Arzela theorem we know that the (/('''s 
admit limit (p in the whole closure of E75. Moreover it is easy to see that cp belong to the set F; 
hence (j) < (j). We necid then to prove only 4> < (p. Let v € F. Similarly to 1), we extend w to u in 
E^ and regularize v to vx in such a way that we have \\v — WA|lioc(2j5) ^ C'A and |V5| < K. Hence 
as before we get \Vvx\ < 1 + CX on E^ and ^ w in C (Ed) as A 0. By simple computation 
we obtain that vx satisfies 

^Avx-\Vvx\^ + l + CX+^Ax>0 in Ed, 
Vx < \x — Qo\ + CX on dTijj, 

where ^a > 0- K we define 

^^A 

Vx 



^l + CX+\Ax 
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by comparison we deduce that 



Choosing d — d'f^ in (31 ) such that 



we see that 

, <4> + CX 

Vi + cx- 

as d'j,^ — ^ oo. Then, letting A — 0, we obtain u < 0; in particular, (j) < 4>. Hence <f> = (f>. □ 



Next we analyze the asymptotic behavior of the solutions of ( 19 1 . From now on in this subsection 
we study separately the two cases < a < | and f < a < tt. Let us consider the first case. 

Proposition 3.6. Suppose that < a < ^ . Let D he a large fixed constant and the solution 
o/(19). Then we have 

^'^ (.t) — > min {di {x) , c?2 (x)} , as d — !> oo, 
uniformly on the compact sets o/S^) H Bd (0), where 



tan^a — 1\ , ,„ / 2 tan a 



di (x) := W X, - , + + U„ - -^—T , (32) 

V V tan a + 1 J \ tan a + 1 ' 



, , . // tan^a — 1\^ , ,,,„ / 2 tan a , 



Remark 3.7. Note that di and d2 are the distance functions, respectively, from the point Qi — 
, , 0, • • • ,0, tan^'^a+i ) ' which is the symmetrical point to Qq with respect to the hypersurface 

tan^a— 1 n n 2 tan a 



ft ^ 

y^tan^a+l'"' '"'tan^a + l^ 

defined by the equation xi tana + Xn — 0, and from the point Q2 — (^\ 



• n • • • n 

tan^a+l'"' ' tan^a + l^ 

which is the symmetrical point to Qq with respect to the hypersurface defined by the equation 
Xi tana — a;„ = 0. So the function (j) (x) is even with respect to the coordinate a;„ and a.e. differ- 
entiable. The problem is that it does not have zero Xn-derivative on {xn = 0}. 



Proof. If 0'* is a solution of ( 21 ) , it is easy to see that 



is a supersolution of ( 19 1 and <j)'^ — sup. 



3.1 



Then must lie in between these two functions. Hence, by Lemma 
the analogous statement for (p'^. The proof of the latter fact is a consequence of Lemma 
the following Lemma [3 



'+sup,gas^ ||x-Qo| + 5 log (U (d {x - Qo))) I 
X — Qo\ + 2 1*^S (d { x — Qo))) I is a subsolution. 

it is sufficient to p rove 
and 



3.5 



□ 



Lemma 3.8. Suppose that < a < ^. If (j){x) is as in (30 1, then 

min{di (x) , ^2 (x)} , 



<i,(x) 



e Bd (0) 



where di and d2 are as in ( 32 I and ( 33 ) 
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Proof. Consider a point x = (xi, • • • , Xn) with Xn > 0. By construction of S^j, the point z G dTju 
which reahzes the infimum will necessarily belong to the set {{xitana + Xn = 0} (1 {xi < 0}}. 
This implies that 

inf Qx - z\ + \z - Qo\) . 

ze{{xi tanQ+x„=0}n{2;i<0}} 

Now we can reason as follows: given x, the level sets of the function z — )• |a; — zj + jz — Qo| 
are the axially symmetric ellipsoids with focal points x and Qq. The smaller is the ellipsoid, 
the smaller is the value of this function; so we are reduced to find the smallest ellipsoid which 
intersects {{xi tana + Xn — 0} D {xi < 0}}. We note that if we fix xi, Xn and vary only x" , the 
corresponding infimum z has the same zi, Zn and different z"; so we can determine zi, z„ in the 
simplest case x" — (0, • • • ,0), and obviously z" = (0, • • • ,0). Then we are reduced to consider the 
minimum problem 



mm 

(zi,z„)e{{a;i tanQ+a:„=0}n{a;i<0}} 



^\/ {xi — zi)'^ + {xn + tan azi)^ + (zi + 1)^ + tan^ azf^ 



Deriving with respect to the variable zi we obtain that at a minimum point 
(xi — zi) + tana (a;„ + tanazi) (zi + 1) + tan^ azi 



{xi — zi)^ + (xn + tan azi)^ \J (zi + 1)^ + tan^ azf 



0, 



which implies 



Zl 



Zn — 



—2 tan axi + (tan^ a — Ij a;„ 
(tan^ a + l) (tanaa;i + Xn — tana) 
2 tan^ axi — tan a (tan^ a — l) .t„ 



(34) 
(35) 



(tan^ a + l) (tanaxi + Xn — tana) 

Now assume that x" 7^ (0, • • • ,0) and xi, Xn are as before. By the previous observation we know 
that the coordinates zi , z„ of the corresponding infimum are given by ( 34 1 and ( 35 1 . So we have 
to determine only z". To do this let us consider the minimum problem 



mm 



(^^ {xi - zif + \x" - z"|2 + (a;„ + tanazi)^ + ^ [z^ + if + |z"|2 + tan^ az^^ 

Again by differentiation we obtain that a minimum point must satisfy 

z"-x" ^ z" _^ 

\J{xi - zif + \x" - z"|2 + (x„ + tanazi)^ ^J (zi + 1)^ + |z"|2 + tan^ azl 

which gives 

„ „ \J {zi + 1)^ + tan^ az^ 

\J (xi — zif + {xn + tanazi)^ + ^ (zi + Vf + tan^ az\ 



(36) 



(37) 



If we plug (34), (35 1 and (37) into (36), we obtain that (x) = di {x). Reasoning in the same way 
for points with Xn < 0, we have 4> (x) = d2 {x). Then we get the conclusion. □ 

Remark 3.9. Note that (j) is a viscosity solution of the Hamilton- J acohi equation |V(/)|2 = 1 in 
In fact, what we have to show is that 
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i) IpP ^ Ij for every x &Ti£i and every p G D^(j) (x), 

a) > 1, for every a; G S/j and every p e (f) (x), 

where D^(j){x) and (f) (x) are respectively the superdifferential and the subdifferential of 4> o,t 
X. Now we can use the description of D'^<p{x) and D~(/){x) given in Theorem 3.4.4 in \CSf : let 
C K" be open and S C ffi*" be compact; let F — F (s, x) be continuous in S x ^ together with its 
partial derivative D^F, and let us define u {x) — niin<,g5 F (s, x); given x G VL, let us set 



M{x)^{s(^S -.u {x) = F (s, a;)} , 
Then, for any a; G 



Y {x) ^ {D,F {s,x) : s G M [x)} . 



and 



D-u{x) 



D+u{x) = co{Y (x)) , 



{p} ifr(x)==p, 

if y (x) is not a singleton. 



(38) 



(39) 



Now we can take fl = 'Ejj, S = {Qi;Q2} and 4'{x) — minjg{i 2} {di (x)}; so 

M (a;) = {g, : (x) = (a;)} , Y [x) = {D,d, (x) : e M (x)} . 



Then, using (38 1 and (39 1, it is easy to see that, if we take x G S^i with a;„ > 0, then D^(j){x) = 
D~(j}{x) — {Dj-di (x)}; in the same way, if Xn < 0, then D'^(j){x) = D~(j){x) = {£'^^2 (x)}. 
So in these two cases properties i), ii) are trivially verified. In the = 0, we have that 

(j) (x) = c?i (x) = d2 (x); then M (x) = {Qi, Q2} arid Y (x) = {D^di (x) , Dxd2 (x)}. Hence, using 
again we obtain D+4>{x) = co{^,^} = ^ 



°{Qi.Q2} 



and D 4>{x) = 0. Then 



we have only to prove property i), since ii) is again trivially verified. To show i) it is sufficient to 
serve that every p £ 0^4) (; 

^2, and that |x — Ql < (x). 



observe that every p £ (x) is of the form p = ^ , where Q belongs to the line joining Qi to 



Let us consider now the case f < a < tt. We have the analogous of the Proposition 3.6 



Proposition 3.10. Suppose that ^ < a < tt. Let D be a large fixed constant and the solution 



of {19). Then we have 



(x) -> $ (x) , as d - 
uniformly on the compact sets o/S/j H Bd (0), where 



$(x) = 



nhn{(ii (x) ,d2 (x)} , 



if tana< ^ii^^^^ 



(l + V^l + ^l) if tana> 



(40) 



(41) 



Proof. We can reason as in the proof of Proposition 3.6 obtaining that it is sufficient to show the 
convergence in (40) for the function 0''. To prove the latter assertion we have to use Lemma 



together with the fact that in the case f < a < tt the function defined in ( 30 ) is equal to that 



3.5 



one defined in (41 1. We can obtain this expression by mixing the arguments used in the proof of 
Lemma 3.8 and those used in Lemma 3.9 in |GMMPl| . □ 
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3.2 Case vr < « < 27r 



In this case we construct the domain S in the following way: we consider the set {a;„ = 0}n{xi < 0} 
and the hypersurface defined by the equation xi tana + Xn — with a;„ < 0. Then we close the 
domain with a smooth surface; the following figure represents a section of the domain in the plane 




We define the scaled domain as in (16) and denote hy To the singularity, which lies on 
{xi = Xn = 0}- As in the previous case, the solution of a Dirichlet problem in S/j will be qualita- 



tively similar to that of ( 15 



We have to study the asymptotic behavior of the solution of the problem 

-^A(/j + (/3 = in Ejj, 



^ = U{d{-^Qa)) on9Ec, 

To do this we consider the function = log<y5, which satisfies 

iA(/.~ |V0|2 + 1 = inSc, 
^^-^log{U{d{--Qo))) ond^D- 



(42) 



Since the asymptotic analysis is very similar to that one made in Subsection 3.1.1 for < a < tt 
we will not repeat the computations. What we obtain is the following result: 

Proposition 3.11. Suppose that tt < a < 2tt. Let D be a large fixed constant and the solution 



of {42). Then we have 

(x) 



dist {x, Qo) = \J (I'l 



as d 



oo. 



uniformly on the compact sets o/S^) H Bd (0). 



3.3 Definition of the approximate solutions 

In order to apply the theory in Subsection |2.H in this subsection we construct a manifold of 



approximate solutions to (10). Since the limit function of the solutions of (19) is not the same 



for different angles a, as we have seen in Subsections |3.1| and |3.2[ we will distinguish the cases. 
We will give the precise construction only for < a < ^; in fact in this case the computations 
are quite different from the flat case a = tt. In the other cases the estimates for the approximate 
solutions are the same (for f < a < tt) or very similar (for tt < a < 27r) to that ones obtained in 
|GMMPl] . Subsection 3.2, and then we will omit the proofs. 
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3.3.1 Case < a < ^ 



Since the function u^ n defined in Subsection 2.2 is an approximate solution of (10) with pure 



Neumann boundary conditions, we need to modify it in the foUowing way. If $ the solution of 



( 19 ), the function 



solves the problem 



= [/(•) ond(ai]c-Qo). 



We can obtain a solution to ( 44 ) looking at the minimum problem 



u 'far oA {\^v\' + v')dy 



(43) 



(44) 



(45) 



From (45) we can derive norm estimate on :^d- In fact, we can take a cut-off function xi ■ 

d (Sd - Qo)^ ]R such that 

Xi (y) = 1 ioT dist{y,d{dJ:D - Qo)) < ^, 
XI (y) = ioryedi^D- Qo) , dist (y, d (SEij - Qo)) > 1, 
|Vxi (y) I < 4 for aU y, 



and then consider the function v (y) — Xi {y)U (y). It is easy to see that 11^11^:^1 



g-d(i+o(i))^ so by ^ we find that 



(d(SD-Qo)) 



< 



(46) 



We can also obtain pointwise estimates on In fact, from Proposition 3.6 we obtain that, as 
d — -t-oo. 



Ed (y) = exp 



\ 



d (tan^ a — l) 
tan^ a + 1 



..II \2 



Vn T 



2c? tan a 
i&v? a + l 



.^o(d)_ 



(47) 

for y G d{V — Qo), where V is any set compactly contained in E^j. Finally, we have pointwise 
estimates for the gradient of S^. Indeed, using the uniform convergence in (20) and reasoning as 



in the proof of Lemmas 3.3 and 3.4 we obtain that ( 29 ) holds true also for Then we can apply 



the arguments in [LNj (see in particular Proposition 1.4, Lemma 1.5 and Lemma B.l) to conclude 
that V^&d — > V0 uniformly as d ^ +00 in any set compactly contained in Sjj on which V0 is 
defined. This convergence implies that, as d — +00, 



VSd {y) = - exp 



d (t&n^ a — 1) 
tan^ a + l 



,,ll\2 



Vn T 



2d tan ( 



tan^ a + 1^ 
V^(^ + go) +0(1)) , (48) 



for y ^ d{V ~ Qo), where V is as before. 
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Now, we want to obtain similar bounds and estimates for ^^^4- and its gradient. Using the 



definition of (y) 



^(1 



Qo) and the fact that also tp depends on d, we have that 



9Sd 
dd 



(y) 



dip /y 
dd \d 



y_ 

d2 



V(y5 



Since p is the solution of (18), we can differentiate (18) obtaining 

in E 



1 A 9ip 



dd 



{x)^WU{d{x-Qo))-{x-Qa) on 91]^ 



(49) 



(50) 



Because of the asymptotic behavior of U at infinity, there exists a positive constant Cd such that 
for d large we have 



1 



C 



D 



u{d{x- go)) < -vu id (x - go)) • [x - go) < CdU {d (x - go)) . 



(51) 



Hence, from ( 18 1 , ( 51 ) , the fact that p > and the maximum principle we obtain that > 
^p in Td- Moreover, as for (19) we can check that the function T'^ := — ^ log<^ satisfies 



dd 



Cd 



T'* (x) 



|yYd|2 

d 



1 



^ = in Sd, 

ilog(-VC/(d(a;-go))-(x-go)) on 



(52) 



Since ^ stays bounded, ^ tends to zero as d 



-oo. Moreover, using again the asymptotic 



behavior of U at infinity, we can say that the boundary datum in ( 52 ) converges in every smooth 
sense (where dTjjj is regular) to |a; — go| as d — )■ +oo. As a consequence, the previous analysis 
adapts to T'' and allows to conclude that still 



and VT'^ 



uniformly as d — )• +oo in any set compactly contained in S^j on which V(f> is defined. 



From (|50|), reasoning as for (|46|), we have that 

+ Qo] 



dp 



< e 



-d(l+o(l)) 



H^d{^D-Qo)) 



3^-v^(I + go) 

9o)- 



(53) 

(54) 
VS, iy) 



On the other hand, from (44) one finds that the function w 
satisfies 

— Aw + ZD = — inc?(I]£) 
d 

To control the boundary value of vu we divide dd{T,D — go) into its intersection with {?/„ = 0} 
and its complement. In the first region we have simply that vu — ^ ' iu)- the second instead 
the estimates in ( 47 ) and ( 48 ) hold true, which shows that the norm of the trace of zu on 



dd{i:D - Qo) is of order e-'^(i+°(i)) 
imply that 



-d\l + 



Mmidi^o-Qo))^^ 



(l+o(l)) 



rf 1+ 



This fact and the latter formula 



Then, from (|54|) and (|55|, we conclude that 



as. 



dd 



(55) 



(56) 



ffi(d(SD-Qo)) 
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Now, using the fact that Lp < Cd \ ^ | and ( 53 ) , together with the Harnack inequahty (which imphes 
IVv^l < Cdif in d (V^ — Qo)) one also finds 



dd 



(y) 



' exp 



yi 



tan a + 1 



+ \y' 



"12 



yn T 



2d tan a 
tan^ a + 1 



1+0 



and 



|V^(y)|<exp 



for d{V ~ Qo) and d — > +oo. 



yi 



2(itan^ a 
tan^ a + 1 



^ , „,„ / 2(itana ^ ^ 

+ y"^+[ynT . 2 ^1 
V tan a + 1 



(57) 



•e°('*\ (58) 



After these prehminaries, we are now in position to introduce our approximate sohitions. Let us 
define two smooth non negative cut-off functions xd : — >■ R, xq ^ satisfying respectively 



XDiy) = l for|y|<f, 
XD{y)^0 for|y|>f, 

(ID 



(59) 



and 



Xo {y) = 1 for y < 0, 

Xq (y) =0 for y > 1, 

Xo is non increasing on M. 



(60) 

Now, using the new coordinates y in Subsection |2.2[ we define 

Ue,Q (y) ■■= Xtio (ey) [{Uq (y) - (y)) (y) + ewQ (y) xo (yi - d)] . (61) 



Following the line of [GMMPlj we prove that the Uc,q's are good approximate solutions to (10) 
for suitable conditions of Q. 



Proposition 3.12. Let /io he the constant appearing in Suhsection \2.2\ Then there exists another 
constant Cn, independent of e, such that, for Co < d < and for Dd < the functions u^^q 
satisfy 



-d E + 



Pi •/2tii 



(62) 



for a fixed C > and for e sufficiently small. 

Proof. Using the coordinates y, we can split u^^q (y) = u^^q (y) + Ue,Q (v), where Ue^q is defined in 
and 



Ue,Q (y) = Xmo (ey) iiXD (y) - l) Uq (y) - XD (y) (y) + e (xo (yi -d)-l) wq (y)] . (63) 
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Then, if we test the gradient of at u^^q on any function v £ Hq (ile), we obtain 
K{'U'e,Q)[v\ = / (V gU^.gV gV + u^,qv) dy ~ ^vdy 

= / gUe,Qy gV + u^,qv) dy - / vF^ Qvdy 
Jvt, Jn^ 

= I'^{u,^q)[v]+Ai+A2, 



(64) 



where 



By Proposition 2.4 and in particular by (14 1 we have that (Wf.g) [v] is of order at most e . Hence 
we only need to estimate Ai and A2 in the last line of ( 64 1 . 



To estimate Ai we divide further Ut^q — ""e.g,! + Ue,Q,2 + u^,q,3, where 

Ue,Q,i (y) = Xmo (ey) (XD (y) - 1) Uq (y) ; m,,q,2 (y) = Xmo (ey) XD (y) Sd (y) ; 

M£,Q,3 (y) = Xmo (ey) e (Xo (yi - c?) - l) (y) . 
Then we write Ai — Ai i + Ai 2 + ^1,3, with 

Ai,i= / {\7gUe^Q^i\7gV + Ue,Q,iv)dy, i = 1,2,3. 
Since xd (y) is identically equal to 1 for |y| < and since xq ivi " d) — \ = Q for j/i < d, from 



([8| and ([Ts]) we get 

To control Ai_2 we write that 

^1,2 = / (Vg?i<,,Q,2VgW + -Ue,Q,2i') rfy = / {g'^^diU^^Q^2djV + u^^Q^2v) dy 

= / {Vue,Q,2yv + Ue,Q,2v)dy + {g'' - S'^) diU^,Q,2djvdy . 

From the condition (c) in Subsection 2.2 we have that — (5*^ | < Ce|y|; then 



(65) 



\A1.2- j {^u,_^Q^2^v + u^^Q^2v)dy\<Ctyj^ |y| |Vt2e,Q,2| dyj ll«llH^(nj ■ 

Since the support of Ue,Q,2 is contained in the set < we obtain from the last formula and 
([46l) that 

(Vw,,Q,2Vi; + ^,,Q,2Z^) dy| < CedDe-'^d+'^d)) • 



1^1 



Now, since satisfies (44), we have 



/ (Vue,Q,2Vt; + -U£,Q,2«) (iy = 
(V (5d (y) (x^o (ey) XD (y) - 1)) V?; + (y) (x^,, (ey) Xd (y) - 1) dy. (66) 
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Since also Dd < c^^i the function (ey) xd {y) — 1 is identically zero in the set < ^} if 
Cq is sufficiently large. Then, using ( 47 ) , ( 48 ) and the Holder inequality, we find that (also for D 
large) 



I / (V (Sd (y) (x,,o (ey) xd (v) - 1)) Vv + (y) {x^, (ey) xd (y) - 1) v) dy\ 



< e 



cLDid / O tan Q(tan Q + . 2d tan 



The last three formulas imply 



1^- 



\ dD , d f 
[T6-+2 V 



SSJ!^ I (l + O 
2 o+lj 



(1)) 



V Ml 



(67) 



D tan n(tan a + 1) . 2d tan a 



From (|65p and the latter formula it follows that 



dDid ! D tan Q(tan , 2d t 



•II^'IIhMo.)- (68) 

It remains to estimate A2. First of all, let us recall that the following inequality holds: 



l<Q -<qI < 



C\Ue,QY' ^\u<l.,q\ for Ue,Q G (O, ^U^.q) , 

C\ue,QY' \ut.Q\ +C\u^,q\^ Otherwise, 



(69) 



for a fixed constant C depending only on p. Moreover, using ([8|) and (13), we can say that there 
exists a small constant ck.u such that 



7 e-^y\ 



1 + 12/1' 



for \y\ < 



1 



We divide next il^ into the two regions 
Bi = < \y\ < min 



d 1 



2' e'^-ff- 



B2 = a\Si. 



For y G Bi we have that x^o ('^2/) = 1, Xr> (y) = 1, Xo (yi - ci) = 1, and hence Ue,Q (y) = -^d (y)- 



By ( |47[) we have also that \u^^q (y) | = |Sd (y) | < e ' x/'-"" "+i ^"^''^ < ifie^Q for y £ Bi. This 
fact, (69 1 and the Holder inequality yield 



\/2d tan c 



Kq - <QMdy < G / K^Qr'\u,,Q\\v\dy < Ce'^-V^^^''"'-'^ Mni^in.) ■ 

On the other hand , in B2 we have that |-u^q| < C {e~i+°'^'^'> +e'~?^^ and that Iw^^q] < 
g-d+o(d). therefore (69 1 and the Holder inequality imply again 

<gll-M2/ < C [(e-^+^f'^) + e-^^) e-+°(^) + e-^^-^'^)] |l.|l,t ^.^^ . 



The last two formulas provide 



\A2\<C 



dp '/2d tan a , ^('^"^ / 



Finally, we obtain the conclusion from (14), (64 1, (68 1 and (70). 



(70) 

□ 
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We have next another estimate for the functional Je, which allows to say that the condition ii) 
in Subsection 2.1 holds true for /. and the manifold of the Uf o's. 



Proposition 3.13. Let fiQ be the constant appearing in Subsection \2.2\ Then there exists another 
constant Cn, independent of e, such that, for Cq < d < and for Dd < -J^^, the functions u^^q 
satisfy 



||/;'Kq)MII < c(^e2 + ee-'^(i+°W)+e"'N''''^^^^^^ ||g|| 

+ C e-^(i+°(i»+e V^^^^iJ' ^'Mll^ll, (71) 



for some fixed C > and for e sufficiently small. In the above formula q represents a vector in 
Hjj (fie) which is tangent to the manifold of the u^,q 's (when Q varies). 



Proof. Since the arguments are quite similar to those in the proof of Proposition |3.12[ we will be 
rather quick. Using the fact that det ((?*■') — 1 and the first line in (64), for any given test function 
V £ Hjj (fij) we can write that 

I'e K,q) M = X! / {9'^c),u,^QdjV + u^^Qv) dy- u^gvdy. 

We want to differentiate next with respect to the parameter Q, taking first a variation qt of the 
point Q for which d stays fixed, namely we take the tangential derivative to the level set of the 
distance d to the interface. Let us notice that in the above formula the dependence on Q is in 
the metric coefficients g^^ and in the function wq appearing in the expression of u^^q (see 61 1. 
Therefore we obtain 



d 



duf 



E 



dg^ 



■diUe^Q^djvdy 



I, J + 



From Remark 



2.3 



(ii) we have that is of order e^\y\. Moreover, computing the expression of 

we obtai^^r%^ = 6x^0 i^y) Xo {Vi - d) gf = o (e^ (l + |y|^) e-l^l), see Subsection 2.2 in 
|GMMP1| . Reasoning as in the proof of Proposition |3.12| we then have 



dQ 



< Ce^ lkll_fi'^(n^) for every v € H\, (fi^) . 



(73) 



On the other hand, when we take a variation qd of Q along the gradient of d, similarly to ( 72 ) 
we get 



dQd 



dUf 



/ 






dQd 






dQ 



diUcQj^djvdy 
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Concerning the derivatives of g^^ with respect to Qd we can argue exactly as for Qt, to find 



Now, computing the derivative of u^^q with respect to Qd is more comphcated than the previous 



case, because ^qq^ has a more involved expression. If we assume that the cut-off function xd iv) 
is defined as xd (f) for some fixed xd, we obtain 



du, 



dQd 



dQd 



eX^oXo (yi - d) 



dwQ 
dQd' 



(75) 



It is easy to see that the last two terms in the right hand side give a contribution to ( 74 ) of order 
at most ee'^^^^"^^)) W^Wh^ (a,) ^^^^ 62grf(i+o(i)) ^^^^ respectively. Concerning the second one, 

we can use the fact that the support of Vxd is contained in the set > j^}, together with 
( 47 ) , ( 48 ) to see that the contribution of this term is at most of order 



( dD.d / P tan Q(taii g + l) , 2d tan a \ / , „ri\\ \ 
\tW + H tan^„ + l + ^,^„,„^J (l + °(l))^^_^(l + o(l)) \ ll^ll^^ 



We can then focus on the first term in the right hand side of (|75|), and consider the quantity 

dd 



/ ( g'^d, ( XmoXd^ ) djv + X^oXi^^w ] dy+p f u^QX^^XD^vdy. (76) 



1.3 + 



Now, using condition (c) in Subsection 2.2 and (56), if we substitute the coefficients g^^ with the 

-d(l+- 



Kronecker symbols wc find a difference of order e e 



^<i(l+o(l)) 



. Next, 



since ^d satisfies — AS^ + S^ = , when we differentiat e w it h re spect to d we get the same equation 

dd 



for so reasoning as for (|66|), (|67|), together with (|57|), (|58|), we find 



It remains to estimate the last term in (76). Usi ng (|56[ ), (57) and the exponential decay of Ue^q 
and reasoning with argument similar to those for (|70p, wc find that it is of order 



g-d(l + 0(l)) +Vtan2„ + 1 



All the above comments yield that 



-d{l+ 



< Cie^ + ee 



o(i)) , -<^/^^^lf5^+7TO](^+°( 



1)) 



+ e 



From ( 73 ) and ( 77 ) we finally obtain the desired conclusion. 



(77) 

□ 
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3.3.2 Case ? < a < tt 



In this subsection we introduce the manifold of approximate solutions in the case f < a < tt. 
Since the construction is substantially the same as in the previous subsection, we will be rather 
sketchy. 



Let us consider the solution of (19 1, and the function defined in (43 1. Reasoning as at 



the beginning of the Subsection 3.3.1 we derive norm estimate for ^d'- 

\\-d\\m(^diSn-Qo)) - ^ 



Moreover, from Proposition |3. 10"| we also obtain pointwise estimates for S^j and its gradient. 

Now, using the cut-off functions (59), (60), we define, in the new coordinates y introduced in 



Subsection 2.2 the functions 

Ue^Q {y) := X^io (ey) [{Uq (y) - (y)) xd (y) + {y) xo (yi - d)] 



Following the line of the Subsection |3.3.1| we prove that the u^^q's are good approximate solutions 
to (10 1 for suitable conditions of Q. Since the computations in the following proposition are the 



same as in Proposition 3.12 and Proposition 3.13 in GMMPT| we will omit the proof. 



Proposition 3.14. Let /Iq be the constant appearing in Subsection \2.2\ Then there exists another 
constant Cn, independent of e, such that, for Co < d < and for Dd < the functions u^^q 
satisfy 



and 



|/"Kq)MII <C^ U' + eexp 



^2 , -rf(l+o(l)) _^g-f^(l+o(l)) 



(78) 



(79) 



for some fixed C > and for e sufficiently small. In (79) q represents a vector in H}^ (rj^) which 
is tangent to the manifold of the u^^q 's (when Q varies). 



3.3.3 Case tt < a < 27r 

In this subsection we introduce the manifold of approximate solutions in the case tt < a < 27r. Also 
in this case we will be very quick, since the construction is the same as in the previous subsections. 



Let us consider the solution of (19 1, and the function ii^; defined in (43). Reasoning as at 



the beginning of the Subsection |3. 3.1 we derive norm estimate for 



P'i|l_H-i(rf(Si,-Qo)) 



< g-'i(l+o(l))_ 



Moreover, from Proposition 3.11 we also obtain pointwise estimates for and its gradient. 



Now, using the cut-off functions (59), (60), we define, in the new coordinates y introduced in 
Subsection |2.2[ the functions 

Ue,Q (y) := X^^o (ey) [{Uq (y) - (y)) xd (y) + ewq {y) xo (yi - d)] . 



Following the line of the Subsection 3.3.1 we obtain that the Ue.q's are good approximate solutions 
to (10) for suitable conditions of Q. Since the computations in the following proposition are very 



similar to those in Proposition 3.12 and Proposition 3.13 in [GMMPl] we will omit the proof. 
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Proposition 3.15. Let /io be the constant appearing in Subsection \2.2\ Then there exists another 
constant Cn, independent of e, such that, for Co < d < and for Dd < the functions u^^q 
satisfy 



and 



{ue,Q) Mil < C (e^ + ecxp-'*(l+°(l» +e-^(i+°(i)) + g-fd+^d))) ||g|| , (81) 



for some fixed C > and for e sufficiently small. In (81) q represents a vector in H]j (fij) which 



is tangent to the manifold of the u^^q 's (when Q varies). 

4 Proof of Theorem 11.11 

To prove our main Theorem we need to derive an expansion in terms of Q and e of tire energy 
of approximate solutions. Then we can apply the abstract theory in Subsection |2 . 1 1 to obtain the 
existence result. 

In the case ^ < a < t: the energy expansions for the approximate solutions u^^q are the same 
as in the case a = tt, see Proposition 4.1 and Proposition 4.2 in [GMMPl] . Then also the definition 
of the critical manifold and the study of the reduced functional are the same. Therefore for the 
proof of Theorem |1.1| in the case | < a < tt we refer the reader to Section 4 in |GMMPlj . 

In the case it < a < 2n, even if the approximate solutions are different from the previous 
case, the energy expansions turn out to be the same. Then also in this case we omit the proof of 
Theorem [O] and refer the reader to Section 4 in jGMMPjj . 

In the case < a < ^ the energy expansions are quite different, so we will give the proof in 
the details. 

From now on we will assume < a < ^ . 

4.1 Energy expansions for the approximate solutions Ue^q 



Here we expand (itj g) in terms of Q and e, where g is the function defined in (61) 
Proposition 4.1. For e — )• and d = d{Q) — > +oo, the following expansion holds 



h {U...Q) =Co- C,eH (eQ) + e-^i+od)) + , V"'" xA^^ J + o (e^) , (82) 

where Cq and Ci are the constants in Proposition 



2.4 



Proof. As in the proof of Proposition 3.12[ let us write Ue,g (y) — u^^q (y) + u^^q (y), see (12) and 



( 63 ) . Then, using the coordinates y introduced in Subsection |2.2[ we find that 

h{Ue,Q) = I^{u^.q)+ {WgU^.QWgile,Q+U^.QUe,Q)dy+l^ (| Vgi2e,g P + W^.g) d?/ 

{\u,,q\p+' - lu^^Q^^) dy (83) 



1 



p+l 



Using condition (c) in Subsection 2.2 we have that 

1/ {V gU^,Q^ gUe,Q + Ue,QUe,Q) dy - I {\/ U^,q\/ lle,Q + Ut,Qile,Q) dy\ 

Jn. JRt 



<Cef |y||V<g||Vii,,g|dy; (84) 
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1/ {\\7gU,^Q\^+ulQ)dy- {\Vu,^Q\^+ulQ)dy\<Ce \y\\Vu,,Q\'dy. (85) 

Jn, JR'i Jr^ 

Concerning (|84|), we can divide the domain of integration into Bd (0) and its complement and use 
([8|, ([13]), (|46f~|47]), ^ to find 



Ce I \y\\W%^Q\\Wu,Q\dy<Ce e-l'^'^''('^^+e-V^ 



For (85), the same estimates yield 



The last two formulas, (83), (84), (85) imply 



le {Ue,Q) = le {Ue,Q) + I {V U^.qW U^^Q + U^.qU^^q) dy + \ [ (| Vm^^q ^ + g) 



+ - 



1 



(86) 



Using the same notation as in the proof of Proposition 3.12 we write Ug^g = Ue,Q,i +We,g,2 + Ue,Q,3- 
Formulas (Isl) and (13) imply 



I / (Vm,,qVz2,,q,i + u,,QW,,Q,i) dy| < Ce-T^r(i+°(i)); 
I / (Vu,,Q Vzi,,Q,3 + Ue,QUe,Q,3) dy\ < Cec-^'^d+'^d)) , 



from which we deduce that 



/ {^U^,Q^U^,Q+Ue,QUe.Q)dy = / {"^Ue,Q^Ue.Q.2 + U^,QUe,Q,2) dy 

Similar estimates also yield 



From a straightforward computation one finds that for any function v 

Vue,Q,2Vw + uI qV = VSd • V (x^o (e-) XDv) + S^x^o (e-) xb"" + V (x^o (eO Xd) (S^Vu - wVSd) 



Applying this relation for v = Me,Q and v = Me,Q,2 respectively, and using ([8|, (13), (46), (47) and 
(lisl) we find that 



+0 e 



(VMe,QVUe,Q,2 + U^,QUe,Q.2) dy ^ I (V (x^o (C') X^M^.q) V5d + Xmo (e-) XDUe.Q^d) 

/ / . I v^dtan a I /-I I „|-1 1 , / dD , d / O t an a ( tan a + 1 ) 2d tan et 

l,'^+ Vtan2a + l J(^ + °(d) + _^ ^-(^Tif + J V ili^^i+I + Vtan2c + 1 , 
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/ (|Vu.,Q,2|' + <Q,2)rf2/= / (\Vix^.o{^-)XD^d)\^ + ix^„ie-)xD^df)dy. 

Using now the fact that, by our construction, the function X/^o (e-) XDUe,Q = Xmo (^') Xd {ut,Q + u^^q) 
vanishes on d (SSj^ — Qo), from (44 1 we obtain 



+ ^ (l^ (^A^o (^•) ^^^-i) I' + (^«' (^•) XuSd)') dy 
(V (Xmo (f) XdUc.q) VS<j + Xmo (e-) XDU£,QSd) dy. 
From ( 86 ) and the last eight formulas we find 

h{Ue,Q) = h{u^^Q) + \ I {\7Ue,Q\/u,_Q+Ue,QUe,Q)dy-\ ^[ (| 

-J^fl+oflll -d- v^"*"" (l+o(l)) 

From (|8|, ([13|, ([14| and ([46| we have that 

J {^Ue^QVue,Q+Ue,QUe,Q)dy = l[ {u^^q) [u^^q] + j \u^^q\p u^^qdy 

and then 

/, K,q) = /, (^2,,q) + 1 [ \u,,Q\Pu,,Qdy + [ (|i2.,Qr+i - K,q\p+') dy 

^ Jo P + i ./o_ 



-d-^M|SJU(l+o(l)) 



ygdtan,, (1 + 0(1)) 
t/tan'^ a + 1 



^ 2g-<i(l+o(l)) 



(87) 



Using a Taylor expansion we can write that 

' ' I o {\u,^q\p\u,^q\ + \u,,q\p+^) otherwise, ^ ' 

As for the estimate of A2 in (70), we divide the domain into the two regions Bi, B2, and deduce 
that 

1 



— T / (KQr'-KQr')dy-- / \u,,Q\Pu,,Qdy 



Therefore using (87) the energy becomes 

le {Ue,Q) = h {Ue,Q) - 



O, 



\/2d tan 1 



I- IV- J.I ^°"^^==f^(l + °(l)) , _ d(p+2) , . 



+ g2g-d(l+o(l)) 
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From (47), the expression of u^ q and the estimates in the same spirit as above one finds that 

\Ue,Q\^Ue,Qdy = - 



and hence from Proposition 2.4 we finahy find 



The conclusion follows from the Schwartz inequality. 



(89) 



□ 



We give also a related result about the computation of the derivative of the energy with respect 
to Q. Again, we will be rather sketchy in the proof since the arguments are quite similar to the 
previous ones. 



Proposition 4.2. For e — >■ and d = d(Q) — )■ +00, the following expansions hold 
^ Ie{u,,Q) = -Cie^WTH{eQ) + o{e^); 

~ n f-<i-4i2jSJLVl + o(l)) , „, 



d 



(90) 
(91) 



where Cq and Ci are the constants in Proposition 



2.4 



Proof. After some elementary calculations, recalling the definition of u^^q in (12), we can write 



I'. K,q) 



d ,_ . f ( _ dii^Q _ du^Q 



where Ue,Q = Ue.Q — We,Q was defined in (63). The first term on the right hand side is estimated in 
Proposition |2.4[ The next two, integrating by parts and using Proposition |2.4[ can be estimated 
in terms of a quantity like 

From the same arguments as in the proof of Proposition |3 . 13| one deduces that the latter integral is 



of order ( 6^2.(1+0(1))^ Vtan2„+i^ ) . To control the first integral in the last line of \92\ 



that 



1+0(1)) {^^ e,-d{i+o{i))-^ 



we can reason as for the estimate of Ai^2 in the proof of Proposition 3.12 to s ee tha t this is of oraer 



From the proof of Proposition 



dQ 



3.13 



one can deduce 



o(e2) +e-3'i(i+°(i)). 



< C (e^ + e ''(i+°(i))) and hence the integral under interest is controlled by 
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Finally, the last term in (92 1 can be estimated using a Taylor expansion as for the term A2 in 
the proof of Proposition 3.12 and up to higher order is given by 

P / Uq~^ (y) ^<i-,Q^Uq (y) ■ qdy, 

where q stands either for the variation of Q in the coordinates y. If q preserves d, the latter integral 
gives a negligible contribution, and we find ( 90 1 . If instead q is directed toward the gradient of d 
the above estimates (and in particular (47l) allow to deduce (91). □ 

4.2 Finite-dimensional reduction and study of the constrained func- 
tional 



We apply now the abstract setting described in Subsection |2.1[ In fact, the following two Lemmas 
hold. 

Lemma 4.3. If Cn is as in the previous section and if we choose 

1 



: Co < d < 



J-<L,Q ■ 



then the properties i), Hi) and iv) in Subsection 2.1 hold true, with 7 — min{l,p— 1}. 

Proof. It is immediate to prove that i) and Hi) hold; in particular, the value of 7 comes from the 
standard properties of Nemitski operators. Property iv) can be easily deduced from the fact that 
the kernel of the linearization of (Hj) in the half space is spanned by . . . , q§^^, as proved in 
|Ohj . and from some localization arguments which can be found in Subsections 4.2, 9.2 and 9.3 of 
[XM] . □ 

Lemma 4.4. For any small positive constant 6, if we take 



u,^Q : (2-(5)|loge| <d< 



eCn 



then also property ii) in Subsection 2.1 holds true, with 



Proof. This lemma simply follows from Propositions |3.12| and |3.13[ 



□ 



As a corollary of the above two lemmas we can apply Proposition |2.1| and Theorem |2.2| so we 
expand next the reduced functional and its gradient on the natural constraint Z^. 



Proposition 4.5. With the choice of Z^^ in Lemma 
have 



4.4 



if is given by Proposition 



2.1 



then 



le {Uc.q) ■■= h {Ue.,Q + We {Ucq)) 



= Co - CieH (eQ) + 6-2^(1+0(1)) ^ 



d[ 1+ 



(l+o(l)) 



(93) 
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d 



d 



-C^e^^TH (eQ) + o (e^) ; 



9 . ^ -d(l+ /St^"° Vl + o(l)) , „, 



(94) 
(95) 



as e — >■ 0, where Cq and Ci are as in Proposition 4.1 and where Qt, Qd o,re as in the proof of 
Proposition |3.13[ 



Proof. By Propositions 2.1 and 3.12 we have that 

Ik. ("e,Q)|| < Ci 11/; (^.,q)|| < C (e^ + ee-'^(i+°(i))) + 



C 



1)) _ d(p+i) 



From the regularity of and Proposition |4.1| we then have 

le {U^,Q + We {Ue,Q)) = /e (Ue,^) + l'^ {u^,q) [We (We.g)] + O (Ue,Q)||^) 

= Co - Ciei/ {eQ) + g-^'^d+'^d)) + e"'('+^teT) + « 



This immediately gives (93), since p > I and (5 is small. 



The remaining two estimates are also rather immediate for p > 2 : in fact in this case property 



Hi) in Subsection 2.1 holds true for 7 = 1, so we also have ||9QU'e|| < Cf (e) by the last statement 
in Proposition 2.1 1. This, together with the Lipschitzianity of I'^ implies that 



d d 

— h iUe,Q) = K iUe,Q + W,) [^qM^.Q + dgWe] = T^h (u.^q) 



dQ 



(96) 



since 7 = 1. The last two estimates then follow from Proposition |4.2[ 

For the case 1 < p < 2, we reason as in the proof of Proposition 4.5 in }GMMPl] to obtain the 
estimates. This concludes the proof. □ 

4.3 Proof of Theorem 11.11 

We use degree theory and the previous expansions. First of all, since Q is non degenerate for H |r, 
we can find a small neighborhood y of Q in F such that ViJ \ry^ on dV and such that in some 
set of coordinates 

deg(VH \r,V,0)^0. 



Then, if S is as in Lemma 4.4 we choose < /3 < |, and consider the set 



y = {(d, g) : d e ((2 - /3) I log e|, (2 + /3) I log e\),eQeV}. 
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Since |r [Q) corresponds to \/tH {eQ) in the scaled domain il^, by using (94) and our choice 
of V we know that, as e — >■ 



Vq^I, (u,,q) = -C^e^VrH (eQ) + o (e^) ^ on -dV. 



(97) 



On the other hand, by (95) we also have 



for d = (2-/3) I log e| 



and 



.^q) = ^Cie^WdH (eQ) + o (e^) , for d = (2 + /?) | log e| 



(98) 



(99) 



Since we are assuming that the gradient of H points toward Sufi near the interface F, V^ff (eQ) 
is negative and therefore the two d-derivatives in the last two formulas have opposite signs. It 



follows from the product formula for the degree and (97)-(99) that 



deg (VI„ y, 0) = - deg (VH |r, V, 0) ^ 0, 

which proves the existence of a critical point for in Y. Since we can choose V and /3 arbitrarily 
small, the solution has the asymptotic behavior required by the theorem, and more precisely by 



Remark 1.2 (6): the uniqueness of the global maximum follows fro the asymptotics of the solution 
and standard elliptic regularity estimates. 



Remark 4.6. To prove also the assertion in Remark 1.2 (a), using (93) in the case of local 
maximum it is easy to construct an open set oj where the maximum o/Ic at the interior is strictly 
larger than the maximum at the boundary. On the other hand, when we have a local minimum, one 
can construct a mountain-pass path connecting the two points parametrized by (7Q, (2 — /?) | loge|) 
and {jQ, (2 + /3) I loge|). Using a suitably truncated pseudo-gradient flow, one can prove that the 
evolution of the path remains inside x ((2 — j3) \ loge|, (2 + /?) | loge|), and still find a critical 
point of Ij . 
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